Quantum phase transitions in bilayer SU(iV) anti-ferromagnets 
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We present a detailed study of the destruction of SU(iV) magnetic order in square lattice bilayer anti- 
ferromagnets using unbiased quantum Monte Carlo numerical simulations and field theoretic techniques. We 
study phase transitions from an SU(iV) Neel state into two distinct quantum disordered "valence-bond" phases: 
a valence-bond liquid (VBL) with no broken symmetries and a lattice-symmetry breaking valence-bond solid 
(VBS) state. For finite inter-layer coupling, the cancellation of Berry phases between the layers has dramatic 
consequences on the two phase transitions: the Neel- VBS transition is first order for all V > 5 accesible in our 
model, whereas the Neel- VBL transition is continuous for N — 2 and first order for N > 4; for N = 3 the 
Neel- VBL transition show no signs of first-order behavior. 



The study of quantum phase transitions is an exciting field 
at the forefront of theoretical condensed matter physics [1|. 
The nature of a particular quantum phase transition is gov- 
erned by properties that affect long distance physics such as 
broken symmetries, topological order, the presence of Berry 
phases and is generally insensitive to microscopic details. 
Quantum magnets provide the richest examples of quantum 
phase transitions because they possess internal symmetries in 
addition to the usual lattice and time reversal symmetries and 
because they often have non-trivial Berry phases in their long 
wavelength descriptions [2|. The most popular internal sym- 
metry group in condensed matter is the SU(iV) group. Initial 
interest was focussed on SU(2) and the case of N > 2 was in- 
troduced purely as a theoretical tool to access the analytically 
solvable N — > oo limit J3][4]]. However, in the ensuing years 
it has come to be recognized that SU(iV) systems with N > 2 
but finite are interesting in their own right, since they serve to 
model a number of physical systems ranging from spin-orbit 
coupled solid-state materials |5| to ultra-cold atoms in optical 
lattice potentials [6|. While the ground states of SU(iV) spin 
models in one-dimensional chains are relatively well under- 
stood [7, 8], two-dimensional phases I9WT1 and their associ- 
ated phase transitions are only poorly understood. 

In this work we address the destruction of the SU(JV) sym- 
metry breaking Neel order in the two-dimensional bilayer sys- 
tem shown in Fig ITJa). In the bilayer geometry the Berry 
phases cancel between the two layers in the continuum limit 
allowing access to the phase transitions of interest without the 
additional complication of quantum interference effects. We 
have studied the properties of the phase transitions from Neel 
order to two different types of paramagnetic states, the va- 
lence bond liquid (VBL) and the valence bond solid (VBS) 
[see Fig.[TJb,c)]. The Neel- VBL transition for N = 2 has 
been studied extensively lfT2l4T5l and is well known to be con- 
tinuous in the 0(3) universality class. Here we address for the 
first time the fate of this transition when N > 2. We find that 
a simple Landau mean-field theory predicts a discontinuous 
Neel- VBL transition for N > 2 and a continuous transitions 
for N — 2. Using unbiased quantum Monte Carlo simula- 
tions we confirm the expectations of the Landau theory, ex- 
cept for N = 3, where we find no evidence for a first-order 
transition. We show that if this transition is continuous, its 



FIG. 1. (a) Bilayer geometry: The white (black) sites are the A(B) 
sub-lattice on which spins transform as the fundamental (conjugate) 
representation of SU(JV). J\ connects nearest neighbors in the plane, 
J 2 connect next nearest neighbors in the plane and J± connect sites 
on different layers. (b,c) shows cartoon product wavefunctions of 
local singlets for the VBL and VBS states. In reality, the ground state 
is a strongly interacting superposition of all valence-bond coverings. 
The ground state nevertheless (b) preserves all symmetries for the 
VBL, but (c) breaks lattice-symmetry (as shown) for the VBS. In this 
paper we provide a detailed study of the Neel- VBL and Neel- VBS 
quantum phase transitions. 



universality class should be identified with a critical point in 
the compact CP 2 model QUEL The Neel- VBS transition 
in the single-layer model has been predicted [ 18 1 and numer- 
ically found to be continuous and in the universality class of 
the non-compact CP^ -1 model for all N |fl"9H2~Tl . We show 
that remarkably the Neel- VBS transition characterized by the 
same broken symmetries becomes first-order in the bilayer ge- 
ometry for all N studied here (our model gives us access to 
N > 5), a striking consequence of the cancellation of Berry 
phases between layers. 

Bilayer Model- Our SU(7V) symmetric model is defined 
with a local Hilbert space of N states on each site of the bi- 
layer square lattice illustrated in Fig. fTTa). We label these 
states as \a) with 1 < a < N. We adopt the representa- 
tion used previously in both analytic EHUI22] and numerical 
ETII23II24I works on bipartite lattices, where the sublattice-A 
states transform under rotations with the fundamental repre- 



sentation of SU(iV) [generated by the N 2 — 1 matrices T a ], 
and the B sublattice states transform with the conjugate of 
this representation. We consider two different SU(iV) invari- 
ant interaction: between sites i and j on the same sub-lattice 
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Using these interactions, we define a 



model SU(iV) symmetric bilayer system as follows, 
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where (ij) denotes nearest neighbors in the square lattice lay- 
ers, ((ij)) denotes next-nearest neighbors in the square lat- 
tice layers and [ij] denotes inter-layer bonds, as illustrated in 
Fig.[TJa). The J\ term by itself gives the familiar single layer 
SU(iV) anti-ferromagnet, which is Neel ordered for N < 4 
and VBS ordered for N > 5. Adding a J 2 term to the J\- 
model favors the Neel state, causing the Neel-VBS transition 
to move to arbitrary large N as J 2 is increased H21H . Finally, 
when the J± term is made large enough it always favors the 
formation of a VBL, by forcing the formation of local singlets 
[see Fig.fTfb)]. The model bilayer anti-ferromagnet, Eq. (fT), 
reduces to the familiar SU(2) bilayer model for N = 2 and 
J 2 = 0. 

Since iJbii satisfies Marshall's sign criteria, it can be simu- 
lated using unbiased quantum Monte Carlo methods on large 
lattices of linear dimension L with 2xLxL sites and at finite- 
temperature T using the stochastic series expansion method 
with loop updates [25-27 1. Neel order is detected by the exis- 
tence of a non-zero spin stiffness p s — T(W 2 ) in the limit of 
L —> 00, where W is the spatial winding number of the world 
lines 1271 . Likewise, long-range order in the correlation func- 
tioniV 2 CV(r,T) = (Po,o + x(0)P r . r+ x(r)) - <P o ,o+x(0)) 2 
signals spontaneous translational symmetry breaking, i.e., the 
onset of VBS order. All the VBS ordering studied in our bi- 
layer system is of the columnar type (at momentum (w, 0)) 
and is in phase between the layers (see Fig. fTTc)). We de- 
fine Oy BS in the usual way as the long distance limit of the 
VBS correlation function. Finally, an absence of both long- 
range Neel and VBS orders indicates the formation of a VBL 
state. Using these tests for the three phases, Neel, VBS and 
VBL, we have computed the T = phase diagram in the 
9± - 92 plane (g± = J±/J\, g 2 = <hl J\) for each N < 10. 
For N < 4, the model Eq. (fTl has only two phases: Neel 
and VBL [Fig. [Tib)]. For N > 5, the model admits in addi- 
tion a VBS phase [Fig.fTJc)]. Phase diagrams for the bilayer 
model, Eq. (fT}, for SU(2), SU(4), SU(6) and SU(8) symme- 
try are shown in Fig. [2] These four cases, contain all the 
types of phase diagrams we have encountered in our study 
with N < 10. We now turn to the main focus of our paper, 
a detailed analysis of the nature of the Neel-VBL and Neel- 
VBS phase transitions that appear in these phase diagrams for 
each N. 

Neel-VBL- First, we analyze the transition between the 
Neel state and the featureless fully symmetric valence bond 
liquid [a cartoon of the VBL state is illustrated in Fig. ITJb)]. 
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FIG. 2. Phase diagram of the model Hui defined in Eq. (fT) for 
SU(2), SU(4), SU(6) and SU(8) symmetry in the plane of g 2 = 
J2/J1 and g± = Jx/Ji. The unfilled symbols are locations of 
first order phase transitions, Neel-VBL (diamonds), Neel-VBS (cir- 
cles) and VBS-VBL (squares). The solid black circles mark con- 
tinuous transitions. For SU(2), the line of Neel-VBL critical points 
shown are in the universality class of the 0(3) non-linear a— model. 
For SU(6) and SU(8) the Neel-VBS transitions shown are in the uni- 
versality class of the non-compact CP^ -1 models (with N — 6,8 
respectively). Solid lines and shaded regions are guides to the eye. 



The Neel-VBL transition in the bilayer model for N = 2 and 
J 2 = has been studied extensively [f2 HT31 . In the special 
case of N — 2 the order parameter describing the SU(2) sym- 
metry breaking can be written as an 0(3) vector. The absence 
of Berry phases in the bilayer geometry then allows for the de- 
scription of the critical point in terms of the well known 0(3) 
non-linear a— model [2|. This simple mapping has no known 
generalization for N > 2. For general N, the simplest de- 
scription of the Neel-VBL phase transition is found by writing 
a Landau theory for the order parameter of the SU(A^) anti- 
ferromagnet. Such a description contains both the Neel and 
VBL phases, since the VBL is featureless and can be thought 
of simply as a phase in which the SU(A^) order parameter is 
quantum disordered. The appropriate order parameter is an 
N x N traceless Hermitian matrix, Q a p, which transforms as 




1.2 1.3 1.4 1.5 

9± 



FIG. 3. Neel-VBL: The spin stiffness p s close to the Neel-VBL 
transition for SU(2), SU(4) and SU(6). The SU(2) transition is con- 
tinuous and in the 0(3) universality class. The quantity p s for SU(4) 
and SU(6) show signs of step-like behavior. Close to the step we find 
double peaked histograms (see Fig. EJ characteristic of a first-order 
transition. The Neel-VBL transition shows such first order behav- 
ior for all N > 4. The parameters used are g 2 = 0.8 for SU(2), 
g 2 = 0.4 for SU(4) and g 2 = 0.6 for SU(6). The legend shows the 
value of L, we have set Ji/3 — L everywhere. 



Q -¥ UQW under SU(iV) rotation. In our model, Eq. (fil 
such a matrix can be constructed microscopically from a local 
operator defined as, Q a p(i) = |a)»(j9|i — 1/N on the A sub- 
lattice and Q a p(i) = |;9)i(a|j — 1/N on the B sub-lattice. We 
can now coarse grain this local operator to obtain the order 
parameter, Q, and write down a Landau theory action, which 
being SU(./V) invariant must consist of traces of powers of Q, 



Sl = oll 



') + /3 L Tr(Q 3 ) + 7iTr(Q 4 ) (2) 



Once the order parameter acquires an expectation value we 
can do an SU(N) rotation to obtain a diagonal form for Q a p = 
™>(fiaifi0i — fiap/N) which is the analog of a "collinear" mag- 
net and the quantity to is the condensate. If we now substitute 
the diagonal form for Q in Sl we can see that generally cubic 



terms in to are present in the action for N > 2. In mean-field 
approach f or N > 2 such terms will render the phase transi- 
tion first order, very much like the first order nematic-isotropic 
transition in liquid crystals |28 ]. When N = 2, it is easy to see 
that Tr(Q 3 ) evaluates to zero and does not give rise to a cu- 
bic to term, making a continuous transition possible. Indeed 
by identifying, n x = (Q 12 + Q 2 i)/2, n y = (Q 12 - Q2i)/2«, 
n z = Qn and including gradient terms in the action we arrive 
at the well known 0(3) <r— model for the n = (n x , n y ,n z ) or- 
der parameter. 

Consistent with the above Landau theory we confirm from 
our numerical simulations (see Figs.[3]l4]and[5]l that the Neel- 
VBL phase transition is continuous for TV = 2 (and in the 
0(3) universality class) and first order for N > 4. The first 
order transitions get progressively weaker as N is lowered. 
Indeed for N = 3 we find no evidence for a discontinuous 
transition up to L < 128 Q2] (see Fig.|5]». If the SU(3) Neel- 
VBL transition is continuous, what is the continuum field the- 
oretic description? Does the field theory admit a critical fixed 
point? The continuum description of the Neel-VBL phase 
transition in our SU(./V) bilayer Hamiltonian is a CP^ -1 field 
theory with a compact U(l) gauge field. In order to make 
this connection, we introduce N complex numbers z a with 
the constraint J2 a \ z a\ 2 = 1 [22 1 and use them to rewrite 
Qa/3 — z a z /3 ~ 5 a p/N . This representation has a well known 
U(l) gauge redundancy which can be made explicit with the 
introduction of a gauge field a M in the long wavelength effec- 
tive action, the famous CP^ -1 model description, 



S = Id xdr 



{(dp - ia^)z a \ 2 + F a pF a p 



(3) 



where F a p = d a ap — dpa a is the EM tensor. Following pre- 
vious work on quantum anti-ferromagnets [29 30 1, it is clear 
that in order for the above field theory to possess the VBL 
state of the bilayer system when J± ^$> J 1; J 2 , the gauge field 
a M must be compact. The Higgs phase with z a condensed cor- 
responds to a phase with SU(A r ) symmetry breaking and we 
identify this phase with the Neel phase. On the other hand, in 
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FIG. 4. Neel-VBL: Hysteresis and double peaked histograms at 
a first order Neel-VBL transition in the SU(6) bilayer. In the main 
frame we show a sample MC history of the binned squared spatial 
winding number, W , which shows clear evidence for metastability 
The inset shows a histogram for the same quantity, with clear dou- 
ble peaked structure. This behavior is found only very close to the 
transition and for sufficiently large volumes, providing unambigu- 
ous evidence for a first order transition. Here shown for L = 32, 
g-2 = 0.6 and g± — 1.36. 




FIG. 5. Neel-VBL: Crossings of the fluctuations of the spatial wind- 
ing number at the Neel-VBL transition for SU(2) and SU(3). In both 
cases up to sizes of L = 128 we see good evidence for a nice cross- 
ing indicating a continuous transition. No evidence for first-order 
behavior was found in these two cases. 
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FIG. 6. Neel-VBS: First order nature of the Neel-VBS transition in 
the two dimensional square lattice bilayer. Both Ovbs an d Ps show 
evidence for step like behavior at the same g±. Close to the jump we 
find the same kind of double peaked behavior in p s that is illustrated 
in Fig. El Here we have shown sample data for N — 8 and g-z = 0.8. 
Similar behavior is found for all N studied here. 



the phase where z a is massive, the photon mode gets confined 
because of the compactness of the gauge field and Polyakov's 
mechanism of monopole proliferation, resulting in a simple 
fully gapped paramagnet, which we identify with the VBL 
phase, Fig. [TJb). Thus the SU(iV) Neel-VBL transition in our 
bilayer can be described in the continuum limit by the Higgs- 
"confined phase" transition in the compact CP^ - 1 theory. Re- 
cent work [16 1 has found that a lattice discretization of the 
compact CP^ -1 field theory has a continuous transition for 
N = 2, 3 and a first order transition for N > 4. Remarkably, 
this is in full agreement with our findings here for the SU(iV) 
bilayer, strengthening the evidence for our identification of a 
continuous transition between Neel and VBL for N = 3. A 
detailed study of critical singularities of the SU(3) Neel-VBL 
fixed point will be presented elsewhere. 

Neel-VBS.- We now turn to the transition between the Neel 
and translational symmetry breaking valence-bond solid state 
[the VBS state is illustrated in Fig. fTTc)]. For a single layer 
the Neel-VBS transition in the model defined by Eq. (fTh was 
found to be continuous ED as predicted by the "deconfined" 
field theoretic arguments [18|. While it is clear that the Neel 
and VBS phases are individually stable to a small but finite 
g±, the interlayer coupling is expected to be strongly relevant 
at the fixed point of decoupled deconfined quantum critical 
points. What is the fate of the Neel-VBS transition in the bi- 
layer geometry? From a theoretical point of view, in the bi- 
layer geometry the cancellation of Berry phases negates the 
quantum interference effects that are crucial to the deconfined 
quantum criticality scenario |[T8ll . In the absence of such ef- 
fects one expects the restoration of the conventional Landau 
paradigm, where the direct transition between two symmetry 
breaking states is necessarily first-order independent of the 
value of N. Indeed as illustrated in Fig. [6] from our QMC 
simulations we find that the Neel-VBS phase transition is al- 



ways first-order in the bilayer geometry. In our model we 
only have access to this transition for N > 5 and in these 
cases we always find a first-order transition. This is a remark- 
able effect since the phase transition in the single layer and 
in the bilayer is in both cases between the same two phases, 
i.e. characterized by exactly the same sets of broken sym- 
metries and in the same spatial dimension. The difference in 
the long-distance physics between the bilayer and single layer, 
much like the Haldane gap in one-dimension, is purely due to 
the presence/cancellation of the Berry phases in the single/bi- 
layer systems. 

In conclusion we have presented a detailed analysis of two 
sets of quantum phase transitions in bilayer SU(iV) spin sys- 
tens: First, we have studied the fate of the popular lfl2l4T5ll 
bilayer S\J(N = 2) Neel-VBL transition for the case N > 2, 
and second, we have studied the fate of the SU( N) Neel-VBS 
deconfined critical point ifTSl I2T1 for a single layer in the bi- 
layer geometry. We have found that the N = 2 continuous 
Neel-VBL phase transition remains continuous for N — 3 (in 
the universality class of the compact CP 2 model 1 16 1), becom- 
ing first order for N > 4, and that the cancellation of Berry 
phases in the bilayer geometry restores Landau's paradigm for 
the Neel-VBS transition, resulting in a first order phase tran- 
sition between two phases with distinct broken symmetries. 
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